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We present results of Monte Carlo study of the monomer-monomer correlation functions, static 
structure factor and asphericity characteristics of a single homopolymer in the coil and globular 
states for three distinct architectures of the chain: ring, open and star. To rationalise the results 
we introduce the dimensionless correlation functions rescaled via the corresponding mean-squared 
distances between monomers. For flexible chains with some architectures these functions exhibit a 
large degree of universality by falling onto a single or several distinct master curves. In the repulsive 
regime, where a stretched exponential times a power law form (de Cloizeaux scaling) can be applied, 
the corresponding exponents S and 9 have been obtained. The exponent S = l/v is found to be 
universal for flexible strongly repulsive coils and in agreement with the theoretical prediction from 
improved higher-order Borel-resummed renormalisation group calculations. The short-distance 
exponents 8 V of an open flexible chain are in a good agreement with the theoretical predictions in 
the strongly repulsive regime also. However, increasing the Kuhn length in relation to the monomer 
size leads to their fast cross-over towards the Gaussian behaviour. Likewise, a strong sensitivity 
of various exponents 9ij on the stiffness of the chain, or on the number of arms in star polymers, 
is observed. The correlation functions in the globular state are found to have a more complicated 
oscillating behaviour and their degree of universality has been reviewed. Average shapes of the 
polymers in terms of the asphericity characteristics, as well as the universal behaviour in the static 
structure factors, have been also investigated. 



PACS numbers: 36.20.-r, 36.20.Ey, 61.25.Hq 



I. INTRODUCTION 

( 2) 

Knowledge of the two-point correlation functions g\j (r) of a polymer chain provides one with profound statistical 
information on the internal structure of the conformation!]] as well as a direct way for computing numerous pair- 
wise observables of interest: mean energy £, mean-squared distances between monomers "D%j, mean-squared radius of 
gyration TZ? g , static structure factor S(p), and many others. While some of these observables permit direct experimental 
measurement and simple evaluation at the level of mean-field theories, the whole range of information contained in 
the pair correlation functions is only accessible in part to the most sophisticated of modern polymer theoriesofj. On 
the other hand, computer simulations based on Monte Carlo Metropolis technique, if performed on a large enough 

(21 

scale for gathering sufficiently large statistics, can yield accurate values of (r,-) on a given spherical mesh n via 
computing histograms of monomer occurrences in spherical layers around other monomers. 

Considerable simulational effort in this areacTQ has been largely inspired by the predictions of the universal scaling 
behaviour-in the field theoretical approaches (e.g_-.hased on the de Gennes' n — > formalism) developed by J. de 
CloizeauxM, B. DuplantieMEl and some othersE3Ej. As simulations have grown larger, the agreement with the 
theoretically computed exponents has improved to a satisfactory leveoj. Many results, however, were only known for 
an open homopolymer chain in the repulsive regime (good solvent), apart from the exact case of the Gaussian chain. 
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In the current work we shall undertake a more systematic study of the pair correlation functions and related 
observables for polymers of varying length in a wide range of the degree of polymerisation N = 100 — 700 units; with 
either open, ring or star (with different numbers of arms / = 3,6,9,12) architectures; with flexible or semi-flexible 
(with stiffness A = 1,5) chain; as well as for the cases of good and poor solvents. 

Generally, the introduction of a more complex chain architecture, stiffness, or attraction, would pose a serious 
challenge for field-theoretic formulations and renormalisation group calculations, and hence only a few scaling results 
are known for the correlation functions in such cases. Moreover, -the SO(n)-</> 4 field theory, which relates to the 
polymer diagrams, was previously only studied up to several loopsEHlj. Recently, however, new results up to seven 
loops have been reported in d = 3 dimensionaU. These, by themselves, would not necessarily yield a better precision 
because of poor asymptotic convergence of the series, but aided by the Borel-resummation renormalisation group 
techniques, the convergence can be dramatically improved and errors in the exponents reduced. Thus, it is important 
to re-examine the agreement of the theory with simulations even in the previously studied case of an open flexible coil, 
as well as to investigate other architectures. Bearing in mind that the theoretical scalings only apply for sufficiently 
long chains, it is also instructive to analyse the data for different chain lengths, and to study dependencies on various 
model parameters generally to address the degree of universality in the scaling laws. 

II. MODEL 

The current coarse-grained homopolymer model is based on the following Hamiltonian (energy functional)0~0 in 
terms of the monomer coordinates, X^: 

H = k ~W E^- X i) J + ^ E A yfe (X l+ X fe -2X,) 2 

+~ £ y(| Xi -x,.|) = i ]T HvdXi-Xjl). (1) 

Here the first term represents the connectivity structure of the polymer with harmonic springs of a given strength 
Kij introduced between any pair of connected monomers (which is denoted by i ~ j). The second term represents 
the bending energy penalty given by the square of the local curvature with a characteristic stiffness A^ between any 
three_consecutively connected monomers (which is denoted by i « j « k) in the form of the Kratky-Harris-Hearst 
termo. Finally, the third term represents pair-wise non-bonded interactions between monomers such as the van der 
Waals forces and so on. In a simple homopolymer model we can adopt the Lennard-Jones of the potential, 

f +oo, r < d 

V(r) = < v /fd\U_(d\*\ _ , , (2) 



M^ 12 -(f) 6 ), r>d 



where there is also a hard core part with the monomer diameter d (below we choose d — £ without any lack of 
generality). „ 

We use the Monte Carlo technique with the standard Metropolis algorithmta, which converges to the Gibbs equi- 
librium ensemble, based upon the implementation described by us in Ref. [If]. Simulations were performed for an 
ensemble of 250 (or 50 — 100 when additional symmetries of the system were present) of independently generated 
initial conditions, each of which has been first equilibrated by typically 20 TV 3 of attempted Monte Carlo steps. To 
ensure proper equilibration, the long- and short-time behaviour of the global observables, such as e.g. the £ and TZg, 
was monitored, as these cease to have a global drift and start exhibiting characteristic fluctuating behaviour around 
well defined mean values upon reaching the equilibrium state for given values of the interaction parameters. After 
reaching the equilibrium, 1000 of statistical measurements for each initial condition have been performed. These 
measurements were separated by typically 40 iV 2 of attempted Monte Carlo steps in between to ensure statistical in- 
dependence of sampling. Thus, overall the averaging statistics was of order Q = 10 5 in all cases. The mean value and 
error of sampling of an observable A are then given by the arithmetic mean (A) = (1/Q) 53? ^7 and by ±y/ (AA) 2 /Q 
respectively. For more details on our Monte Carlo technique we refer to the previous papers Refs 

III. DEFINITIONS 

The intra-chain correlation function of a pair of monomers i and j is defined as, 



id by ± 

HE 



2 



,(2), ^ / — „ ^\ I 



<^(r) ee ^(X, - X, - r)^ = ^fflX, - X,| - r)j. (3) 

The second equation establishes that it is a function of radius r = |r| only due to spatial isotropy (SO(3) rotational 
symmetry). We may note that this function should, strictly speaking, be named distribution function, but since 

f2l (2) 

9ij ( r ) ~* ^ when r — > oo because of the chain connectivity, we apply the term 'correlation function' to g\^ (r) itself 

rather than to the quantity g$ (r)/(g^) 2 — 1, which would vanish as r — ► oo in the case of simple liquids. The 

function is normalised to unity via: J d 3 r g?) (r) = 1 . Note that the correlation functions exactly satisfy the excluded 
volume condition, 

g ( i f{r) = forr<d, (4) 

due to the choice of the hard-core part in the non-bonded potential Eq. (H). 

The mean energy and the static structure factor (SSF) then can be computed as follows, 

£ = = " £ |d 3 ri^(|r|)«4 2) (r), (5) 



i 1 r°° 

S(p) = jjY,~9 i2) (\P\)> 5 (2) (p) = <exp( l p(X,-X J )) = ^ J a r 2 dr 



sin(pr) ( 2 ) 



pr ^ 



9\?(r), (6) 



where tilde indicates the 3-d Fourier transform. The mean-squared distance between monomers i and j is, 



T>ii = (Ai) - ^(X,-X,)^ = I d 3 r|r| 2 ffl ( f(r), 



(7) 



which we defined here without the traditional factor of 1 /3 as compared to previous paper Ref . [l'J. Then the mean- 
squared radius of gyration is simply, 

ij, i^J 

Finally, for analysis of the polymer shape we can define the non-averaged shape tensor0~0, 

Q^OQ = 2^2 E W - *?)(*f " (9) 

which is related to the inertia tensor 7 Q/3 = (5 Q/3 tr Q — Q Q/3 in the case when all monomers have equal masses. Let 
us also introduce the traceless tensor Q a ^ = Q a @ — ^S a ^tiQ and the 'mean' value of the eigenvalues q = (l/3)trQ, 
where the trace of Q according to Eqs. (|]]|) coincides with the non-averaged squared radius of gyration, namely: 
tr Q = So=i 2 3 1^ = Rg- One can study the mean ratios of the eigenvaluesEH, 

* w -<£>-<£>■ .£/ w ^ 

In addition, we can define the following two asphericity characteristics, &0, 

As = {As) = = ( J E 3 iqia) e q? ) ' ( s P here ) ^ ^ 1 Mlmear), (11) 

/27detQ\ / ria = i,2,3(g (a) -g) \ l,, lfWC ,.v noA 

Ss = (Ss) = ( -T^QjF > = ( ^3 )> -4 (oblate) < 5 3 < 2 (prolate). (12) 

Alternatively, we can consider modified quantities ,0011 

A {fAs) S {fSs) (13) 

A * = TFT' 3 = ~¥V' { ' 

which are more amenable to analytical treatments, although less sensitive on the shape of conformations. 
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IV. CORRELATION FUNCTIONS 



A. Flexible ring polymers 



Ring polymers are the simplest objects of study since they possess exact translational invariance along the chain, 
so that all single-point observables are constants and all pair-wise observables are functions of the separation in the 
chain indices n = \i — j\ rather than of both indices0, i.e. T>ij — V$ n , g^\r) — (r). 

In Fig. [l] we exhibit the mean-squared distances T> n of flexible ring homopolymers vs the separation n. There is 
an additional symmetry property of rings such that T> n is also symmetric with respect to the middle (i.e. maximal 
separation along the chain), namely T> n — T>N- n - This function has a characteristic bell-shape (which would become 
somewhat deformed in case of semi- flexible polymers), and this shape agrees very well with that from the Gaussian 
Self-Consistent (GSC) method (see e.g. Fig. 2a in Ref. [l6]). This function also appears to possess an approximate 
scaling transformation property, so that curves with different N can be superimposed into each other by the trans- 
formations: n — n/N, T>n — T> n /N 2vF , where vp ss 3/5 is the value of the Flory swelling exponent. Such rescaled 
quantities T>n are shown in Fig. || for both the coil (in a good solvent) and the globule (in a poor solvent). For 
the coil, T>n coincide with each other very closely everywhere except around the middle of the range |0.5 — h\ < 0.2, 
where such agreement is less accurate. However, without rescaling, T> n is strongly dependent on both n and N. The 
dependence on N is relatively weak only for n corresponding to a few polymer links, where connectivity is directly 
manifested. 

("2") 

Similarly, without rescaling, the intra-chain pair correlation functions g n (r) are strongly dependent on both n and 
N. However, if we introduce the rescaled pair correlation function in terms of the dimensionless variables, 

a (2) /~n _ ^3/2 (2) / n. a _ T . 

9l j '(r) = V l < gy>(r), r =tfj2> ( 14 ) 

ij 

(2) 

these would permit much more straightforward comparison with each other. Note that g\j (f) satisfies the following 
two normalisation conditions: 

i rOO i 

drr 2 gg\r)= / drf i g ( i f(f) = —. (15) 



Thus, in Fig. we present §n (f) for n corresponding to the half-ring chain separation for different values of 
the degree of polymerisation N. It is important to emphasise that the dashed curve here corresponds to a different 
connectivity constant and this will be discussed later on. The two lower curves coincide with each other remarkably 
well, except for very small values of f, which are affected by the hard-core part of volume interactions in the area 
of direct steric repulsion r ~ d (see also Eq. (^) . The same is true of gff (f) curves corresponding to different chain 
separations n, except very small values of n < ~ 5 for which T> n is of order of several d 2 units, i.e. where T> n d 2 > 1. 
We may remark that at such very small n any particular polymer model (e.g. bead-and-anharmonic-springs one 
or a freely-jointed one) would in any case exhibit its model-specific features, and hence would no longer produce a 
universal behaviour. 

One can see that the function g^ 2 \f) is non-monotonic: it increases for small f, reaches a maximum at around 
f ~ 1/2 and then decreases in a Gaussiaprlike manner, being almost zero for r > 2. The short-distance behaviour, 
which is often called the correlation hola3, is an effect of the excluded volume interactions coming from repulsions 
mediated through higher than direct binary contacts described by Eq. (^). 

From this interpretation one would expect that the correlation hole size should decrease upon decreasing the 
excluded volume size d relative to the Kuhn length i/y/R. This is indeed the case as can be seen from the dashed 
curve in Fig. |^, which has been obtained for a polymer with a twice weaker spring constant k. The correlation hole 
size decreases, and the value of g^ there correspondingly increases to keep the normalisation condition and hence the 
function reaches the maximum earlier. The effect is the opposite at distances f ~ 1, but much weaker, and the tail at 
larger f is practically unaffected by this change of connectivity parameter. We would like to emphasise that while the 
change in the correlation hole with the Kuhn length is quite dramatic, it is by no means unexpected. Evidently, as 
much weaker springs are considered (with d fixed) , we should approach the Gaussian law for the correlation function, 

gn^ oc exp(— 3 f 2 /2), as well as for the mean-squared distances, and the effect of cross-over behaviour is indeed so 
clearly seen in Fig. |[ 
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B. Globule of a ring chain 



After collapse of the coil, the mean-squared distances T> n of the globule at the van der Waals attraction Vb = 6 fe^T 
are shown in Fig. |J vs the chain separation n for different chain lengths N. These functions increase rapidly at small 
n, quickly reaching a stable plateau scaling as V n ~ N 2 / 3 at n ~ n* ~ N 2 ^ 3 . The plateau is related to the nearly 
constant density within the globule. The shapes of V n in Fig. || agree with those from the GSC method (see solid 
line in Fig. 2a in Ref. |l6| and curve marked LLG in Fig. 2 in Ref . E9) . 

Rescaled curves T>n = N~ 2 / 3 T>n vs h — n/N in Fig. || lie fairly close to each other. With increasing N the plateau 
height in T>n slightly decreases, reaching its asymptotic value as the cross-over scale vanishes: h* ~ TV -1 / 3 — > 0, so 
that the normalisation Y^n<i = N~ 2 / 3 remains constant. 

In Fig. H we present the rescaled correlation functions fjn\f) of the globule for several different values of n. These, 
for small f, have several peaks of increasing width and decreasing height located at integer multiples of the excluded 
volume diameter r = dn, which is similar to the behaviour of correlation functions in simple liquidsc3. For large f this 
behaviour changes into a monotonically decreasing tail. All functions gn (?) for chain separations \N/2 — n\ > n*, 

( 2) 

where n* is the same cross-over index as in V n , lie very close to each other and to the half-ring function 9 n / 2 (t)- 
Thus, as long as the chain separations n are not directly affected by the chain connectivity the idea of a single master 
function in terms of the dimensionless variables is still valid for the globule, although the function is much more 

(2) (2) 

complicated than for the coil. This is consistent with the fact that peaks in g„ (r) occur at r = dn, so that for gk '(f) 
to fall onto a single function T> n has to become n- independent, which indeed happens beyond n*. 

(2) 

Importantly, the curves gn do not fall onto a single master curve for different chain lengths N here. Indeed, 

-s (2) 

the number of peaks in g n can be estimated as proportional to the linear size of the globule, which is of order 
oc N 1 / 3 . Thus, for a very large N we shall have a function containing a sharp and narrow first peak followed by 
many oscillations decreasing in amplitude around a constant level (related to the constancy of the density inside the 
globule) and ended by a monotonically decreasing tail. 



C. Open chain 



Rescaled correlation functions of an open chain are presented in Fig. ]a. The ring symmetry no longer applies here, so 

(2) f2l (21 

that g\j (f) now explicitly depends on both chain indices. We note that the reflection symmetry g N+ i_, t N+i-j = 9ij 
is the only exact symmetry of an open chain. Thus, in obtaining the correlation functions, we may average over two 
possibilities for the end-internal and internal-internal correlations. As for the end-end correlation function only single 
possibility exists hence leading to \/2 times worse statistics in this case. Since the averaging statistics is the worst for 
an open chain as compared to other architectures considered in this paper, we present the correlation functions here 
via data points and smooth curves approximating them (this will be discussed in more detail in Sec. [v|). 

Clearly, the end monomers have a special role due to entropic reasons, but correlations of internal monomers on 
large separations within a very long chain should behave in a similar manner to those of an equivalent ring. Indeed, 
the function g^j A 3^/4^) (denoted by dotted line in Fig. ^) practically coincides with that of a ring (solid line in Fig. 

I)- 

However, the functions involving the end monomers have distinct shapes. Short-dashed curve of N j 2 {^) nas a 

higher peak and a somewhat smaller correlation hole area than the middle-middle correlation function; the trend con- 

(2) 

tinues for the long-dashed curve of 3^/4(^)5 w hich shows a weak cross-over behaviour; and the end-end correlation 

function g\ N (f) (solid line) has the tallest peak and the narrowest correlation hole area. Interestingly for f > 1.5 all 
functions still have practically identical monotonically decreasing tails. 



D. Semi— flexible ring 



Now then, it would be interesting to understand the effect of chain stiffness determined by the Xijk constants in 
Eq. (0). Thus, let us consider a semi-flexible ring polymer of N = 200 units with A^-fc = 1 for each triplet of 

consequently connected monomers and zero otherwisetS In Fig. ^ we present g^\f) functions for different values of 

chain separation n. As opposed to similar functions of a flexible ring (A = 0) in Fig. |3[ in this case gn\r) no longer 
fall onto a single master curve. 
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A striking feature for relatively small n is the expansion of the correlation hole area at small f. The latter is now 
determined not only by the excluded volume interactions but also by the tendency of the chain to remain locally 
straight via the stiffness term. The lowest curve for n = 10 in Fig. [7] shows that this effect is still present for n 
considerably larger than the persistence length of the chain, which is of order A. Upon increasing n the correlation 
hole shrinks, producing a steeper growth at small f, and for large enough n the hole becomes even smaller than that 
of the flexible ring. This may be understood as follows. Due to the stiffness effect neighbouring monomers along 
the chain are being pushed away from the current monomer, thus providing an easier access for contacts with it for 
monomers which are more distant along the chain. The mismatch of g n for different n is present also at intermediate 
values of r, where a complicated cross-over occurs, and, unusually, even the tails of the distributions do not match 
perfectly here. Thus, stiffness completely removes the universality of rescaled correlation functions which was present 
in the flexible-.case for topologically equivalent monomers. Moreover, manifestations of stiffness are also strongly 
iV-dependentiij and the consideration of the limit N — > 00 is highly non-trivial as it would require renormalisation of 
A, which might significantly change the physical persistent length of the chain. 



E. Star polymers 

Now let us turn our attention to the case of flexible homopolymer stars with arm length (N — 1)/ f = 50. We refer 
the reader for notations, further details and previous results on stars to our preceding work in Ref. [l^. 

Clearly, an open chain may be viewed as a star with / = 2 arms, so one would like to understand the effect of 
increasing the arms number, which will be studied in Sec. but we shall first analyse different correlation functions 
in a star with sufficiently large /. Thus, in Fig. ||we present the functions for different ij for a star with / = 12 
arms in the coil state. 

Star architecture allows us to use the arms symmetry to increase averaging statistics. Therefore, it is sufficient to 
consider correlations of three types: core-arm, intra-arm and inter-arm. As can be seen from Fig. @, the shapes of 

functions g^-' corresponding to these three types are very different indeed. 

The core monomer plays a very special role, so the solid (core-half arm) and long-dashed (core-end monomer) 

curves in Fig. M show pronounced correlation holes at small f due to stronger steric repulsion from the core monomer 

(2) 

as / increases. Interestingly, g^ grows faster at small f in the core-end case than in the core-half case. This may be 
explained by entropic reasons as it is easier for an end monomer to form a single long loop and return towards the 
core. 

The intra-arm correlation function drawn as short-dashed curve in Fig. pi has a shape and meaning similar to that 
of the corresponding function in an open chain (short-dashed curve in Fig. ^J) since the core monomer and the rest of 
the star are typically far from the part of the arm in question. Analogously, the inter-arm end-end correlation function 
(dotted curve in Fig. ||) has a similar shape and meaning to that of the end-end function in an open chain (solid curve 
in Fig. U), although the correlation hole is smaller in the latter case. This is because there are more possibilities for 
any two ends in a star to encounter each other than to do so for the only two ends of an open chain. This effect will 
be numerically investigated for different / in Sec. |v| Thus, correlation functions of a star are essentially ij-dependent 
even within monomers of identical topological types (i.e. inner, end or core). 



V. SCALING RELATIONS 
A. Open and ring chains 

In this section we shall quantify the behaviour of rescaled correlation functions in a good solvent. According to 

Refs. |l]jL3| these, at least in the well known case of an open flexible chain, can be described as a power law times a 

— (2) 
stretched exponential, so that we can fit (f) there via the following function, 

gg\r) = A ij r e » expi-Bijf*). (16) 



For brevity we shall suppress the ij indices bellow. Due to the two normalisation conditions in Eq. (15) constants A 
and B can be immediately calculated and expressed via 9 and 5: 

r((5 + 6)/5) SB^/ S 
r((3 + 8)/Sy/ 2 ' 4tiT((3 + 0)/*)" 1 ' 



G 



The exponents S and which in the case of end-end correlations of an open chain is denoted as #0j can be expressed 
via, 

1 v- 1 

s = - — , e = l — , (is) 

1 — v v 

where v has the meaning-of the inverse fractal dimension of the system and 7 is related to the number of different 
polymer conformations™^. 

We should, strictly speaking, comment at this stage that for small distances the #0 exponent should, in fact, be 
replaced by another value ae = (1 — 7 + vd — d/2)/(l — v) = 0.249 ± 0.011, but because ae is numerically very close to 
9q in d = 3 dimensions, one .can use Eqs. (Jlfjjl^) in a single form for all distances. Such formula is often called the 
des Cloizeaux scaling forrrotij. This expression also works very well for other types of correlations and we shall fit 
our data via Eq. (Hq) below. 

The exponents S and v have been derived in Refs. |27]j28|. The theoretical values that we shall quote in this paper 
have been updated by a higher order calculations in Rcf. |14|, 

7 = 1.1596 ±0.0020, v = 0.5882 ±0.0011. (19) 

The exponents V , v — 0,1,2 related to the probabilities of contacts of the end-end, end internal and internal- 
internal monomers in an open chain respectively, have been calculated by des Cloizeaux in Ref. |^ by the renormalisation 
group technique based on e = 4 — d expansion for the Euclidean field theory, as well as by Duplantier in Ref. |TI] and 
Oono et al in Ref. |l2| based on the Edwards' model. 

In Tab. | we present the exponents, based on fitting of our Monte Carlo data, as well as the best theoretical 
values known to us (the latter and the corresponding Monte Carlo values under comparison are_aet in bold face). 
Fitting has been done via the the nonlinear least-squares (NLLS) Marquardt-Levenberg methodcj by means of the 
fit function in the gnuplot software ver. 3.7.0. Fit reports parameter error estimates which are obtained from the 
variance-covariance matrix after the final iteration. By convention, these estimates are called 'standard errors' and 
they have been reported in the tables. 

Two types of fittings have been independently carried out. In one case, both 5 and 9 have been fitted, while in the 
other case, S was fixed at the theoretical value and hence a more accurate estimate has been achieved. The latter 
exponent will be denoted as 9s fix ■ 

One should bear in mind the reservations that universal exponents values are only asymptotically approached as 
TV increases, and that assuming the same power law for both small and large distances is only approximate. Yet, we 
find our results in a good agreement with the best of theoretical estimates. 

First of all, the exponent 8 in cases when fitting errors are the smallest seems quite close to Stheor- The agreement 
is particularly good in case 1', where 9g fi X is also very close to theoretical value 9\. This may be in part due to 
the fact that here the two internal monomers are well separated within the chain, reducing finite-size behaviour. It 
is interesting that S for the ring is also close to the current exponent, which indicates that as N increases the same 
universal exponent value should be achieved in this case also. We should comment that a somewhat longer ring 
with N = 300 is equivalent to only a N — 150 open chain though. As for the end-end case 0, it is somewhat more 
intrinsically noisy, since ends have more entropic freedom, so that the agreement with the much more accurately 
known 9q is only satisfactory. 

In the internal-internal case 2 our fitting error is rather minimal, thus our value of 9s fix here is also very consistent 
with the ten times more accurate result for the ring. Such improved accuracy is achieved for the ring because here 
one can average over all pairs of monomers separated by n = \i — j\ links due to the kinematic symmetry described 
above. Thus, the error is reduced by an order of vT00, which is also seen in the data. It is only natural that any 
two well separated internal monomers in an infinitely long open chain should behave as those of an infinitely long 
equivalent ring since an extra link between the two ends should not matter. However, the theoretical estimate for 92 
is the least accurate of all, so that higher order calculations akin to those done for 9o recently may be required here 
too. 

Moreover, here we can see a cross-over effect discussed in Ref. || at work. By taking one monomer at the end, and 
by moving another from 1/4-chain (case 1") to 1/2-chain (case 1) and to 3/4-chain (case 1'), the exponent 9s fix 
would steadily decrease, and as we reach the other end (case 0), the smallest value 9o would be reached. The tilde in 
front of the theoretical exponents in Tab. | indicates that those exponents are in the cross-over area. 

B. Stars 

Results from fitting of various correlations functions for stars with different / in the manner described above are 
collected in Tab. O. In this case, no theoretical exponents are known to us as universality is not a feature of this 
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system. Yet, it is instructive to attempt the fitting even though we realise the finite-size and cross-over nature of the 
resulting numbers. This should not be very off-putting though as actually synthesized stars are indeed objects of a 
few arms of very finite length at present. 

First of all, we find again that the stretching exponent 8 agrees within fitting errors with the theoretical value for 
open and ring chains Stheor- Thus, we can concentrate on the discussion of a more accurate exponent, 8$ fi X . Clearly, 
the core monomer has a very special role, and as the number of arms grows, the steric repulsion around it becomes 
really significant. This leads to a rapid expansion of the correlation hole as / changes from 3 to 12, so that values of 
@8 fix grow rapidly with / for the correlations functions of the core with the middle and ends of arms (two first rows) . 
As for the intra-arm correlations generally, these behave more and more as in an open chain as long as the monomers 
in question are away from the core. So the middle-end intra-arm exponents (third row) 8$ fi X are fairly independent 
of / and are close in numerical values to the corresponding exponent 8i of an open chain. 

Now let us discuss the inter-arm exponents. The middle-middle (fourth row) and end-end (last row) inter-arm 
exponents 8s fix start from values close to those of the corresponding exponents of an open chain 82 and 8q respectively 
at / = 3, but both of them decrease significantly as number of arms / increases. Finally, the end-middle inter-arm 
exponent 8s fix (fifth row) starts from a value fairly close to 8\ and decreases only slowly with increasing /. 



VI. ASPHERICITY CHARACTERISTICS 



In this section we would be interested to investigate the average shape of conformations, for which we have studied 
the spherically averaged pair correlation functions above. 

The shape tensor in Eq. (^|) is merely a slight generalisation of its 3-d trace 7Z. 2 , in Eq. (||). However, since the space 
rotational isotropy is enforced on average, one way of obtaining a non-trivial spatial information is to diagonalise the 
matrix Q a @ for each member of the statistical ensemble and then to average some useful observables constructed out 
of its eigenvalues. Eq. (|l0|) introduces the averages of normalised sorted eigenvalues (so that A^ 1 -* > A'- 2 -' > A' 3 '), 
whereas Eqs. ( |Tl| ) and ( |12| ) introduce the asphericity index A3 and the ellipsoid index S3. The latter allows one to 
further distinguish a prolate ellipsoid (rod-like) from an oblate one (pancake-like). Finally, variables A3 and i^pin 
Eqs. ( [l3"l ) introduce related fractions of averages, which although less informative are easier to compute analytically!! 3 ! 

In Tab. Ill we have collected these asphericity characteristics for the systems under study. First of all, we may 
observe that the 'ellipsoid aspect ratios' A^ a ' reflect the shape fairly well, although the asphericity index ^3 and the 
ellipsoid index S3 express it in an even more clear way. As for the observables A3 and S3, they follow the trends of 
A3 and S3 respectively, although in a somewhat less sensitive manner. 

Now then, for the ring, open coil and the globule we see that the values depend little on N, if the system size is large 
enough, approaching their universal limits. The coil of a ring (first row) has a smaller asphericity and prolateness 
(rod-like shape) than the corresponding values for an open chain (third row) . We may note also that our values for 
an open coil are fairly close to those reported in Refs. ^3|,|l3[ The semi-stiff ring (fourth row) becomes, clearly, more 
aspherical and prolate as stiffness A increases at fixed N. If, however, A is kept fixed, and N increases both A3 and 
53 decrease approaching the flexible chain limit since the ratio X/N presents the relevant degree of flexibility of the 
polymer. 

The globule of a ring (second row) is nearly perfectly spherical and only very slightly prolate, as should be expected, 
and these asymmetries diminish even further as N increases. However, at the same time, the non-equality of A^ a ^ is 
more noticeable here, which shows that they provide less informative asphericity observables. 

Finally, stars (fifth row) starting from a fairly non-spherical and prolate shape at / = 3 become significantly more 
spherical and less prolate as the number of arms / increases. This, of course, is consistent with the intuitive idea that 
a star with a very large number of relatively short arms behaves as a sphere due to strong steric repulsions. 



VII. STATIC STRUCTURE FACTORS 



The static structure factor (SSF) is of interest due to its relation to the light and neutron scattering techniques and 
it is standard to introduce it in the rescaled Kratky form via, 

S(p) = ^p L S(p), P = P^ g , (20) 

which has an advantage of being less sensitive on N. _ 

The factor p 2 allows one to emphasise the large p behaviour more clearly. Indeed, it is well knownEj that S(p) cx 
p~ x / v with v = 3/5 and 1/2 for the repulsive and ideal coils respectively. Thus, S(p) would be increasing, reaching 
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a constant asymptote, or decreasing for the repulsive coil, ideal coil and globule respectively. This is indeed the case 
in Fig. [)[ where we depict S(p) vs p for various systems. Since curves for the coils (or the globule) for different ./V 
coincided with each other so well due to a nearly perfect universality, we only have drawn one particular size in each 
case. i — | 
Clearly, for small p we have the identical behaviourcJ, 

S(p)^f(l-p 2 /3 + ...), (21) 

in terms of the rescaled variables in all cases. The SSF of the globule has a characteristic oscillating behaviour, which 
shows its dense structure. It is instructive to compare this with the SSF of a solid sphered of a radius R s , 



S( P ) = (j^jE [sm(RsP) ~ (R.P) cos(R sP )j 1 . (22) 

By requiring that the cumulant law Eq. ( pl| ) is universal in terms of the rescaled variables one determines the radius 
to be R s — J (5/3)7£g. One can see in Fig. || that the SSF of the globule is in a perfect agreement with that of 

the solid sphere up to the peak point. At larger p both functions have rather similar shape, however the SSF of the 
globule has somewhat smaller amplitude of oscillations, which reflects fluctuations in shape and size of the globule. 
As additional analysis shows, with increasing N the minima of the rescaled SSF of the globule decrease in value, so 
that S(p) asymptotically approaches that of the solid sphere. 

The SSF S(p) of an open coil has a power-law increasing behaviour, which is preserved for a ring coil as well, 



although the latter has a slight peak, showing a more spherical structure of a ring in accord with Sec. VI. However, 
when we introduce even a moderate stiffness for a ring coil, the plot of S(p) starts to significantly deviate from that 
of a flexible ring, and also becomes quite N dependent. This results from the loss of master function unive rsality in 



the rescaled correlation functions in terms of N, which we have seen for a semi-flexible coil in subsec. IV D . 

Now, if we turn our attention to the structure factors of stars with different numbers of arms in Fig. |L0[, we can 
see how the starting case of / = 3, which is quite reminiscent of the curve of an open chain, transforms into a shape 



corresponding to more spherical objects as / increases (as we have also seen in Sec. VI). Again, the /-dependence 



here is quite dramatic, although S(p) still increases for large p as the stars are in an extended state. 



VIII. CONCLUSION 



In this paper we have studied the correlation functions and related observables for a single homopolymer chain with 
several macromolecular architectures. 

A ring case is the most symmetric one and here we have discovered the universality of the dimensionless correlation 
functions gn^ for all chain index separations beyond small scale w) ~ 5. Here the universality also holds for different 
chain lengths N. 

The universality appears to be valid in terms of n for the globular state as well, although the scale index is of order 
n* oc iV 2 / 3 there, and the correlation functions have a complicated oscillating behaviour. Interestingly, the original 
function g\ (r) is also n-universal for the globule, which is consistent with the behaviour of the pair correlation 

(2) 

function akin to that of simple liquids, so that the peaks in g n (r) occur at integer multiples of the hard core 
diameter r = md, requiring for T> n to become constant beyond n* . 

For the coil the functions g„ can be accurately fitted via a power law times a stretched exponential — the des 
Cloizeaux scaling function Eq. (|l6|). The stretching exponent 8, related to the fractal dimension of the chain, as 
well as the power exponent 9, related to the correlation hole effect, are found to be close to the results of theoretical 
calculations in the special case when the Kuhn length is chosen equal to the hard sphere diameter: £/(f/nd) — 1. This 
strong repulsion regime is a special one because it most closely corresponds to the continuous model of the polymer 
chain. 

However, we have found that decreasing k results in a significant change in the power law exponent 9, even though 
the long distance exponent is less affected by that. This may be easily understood as a cross over to the Gaussian 
chain behaviour, namely 9 — > and 8 — > 2 as £/(y/Kd) — ► oo. Strong sensitivity of the correlation hole on the model 
parameters also demonstrates that care should be taken when comparing universal renormalisation group predictions 
with results from simulations of concrete systems. 

Clearly, the ends play a special role for an open chain. Therefore, one has to distinguish the end-end, end-internal 
and internal-internal correlations with the corresponding exponents 9 V respectively. The exponent 8 is found to be 
universal for strongly-repulsive flexible coils of any architecture and numerically very close to that of the ring chain 
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and to the theoretical prediction. Fitting '(f) via the des Cloizeaux expression in the strong repulsion regime 
(£/(y/K,d) = 1) has produced values fairly close to the results of theoretical calculations. The effect of cross-over, 
when the internal monomer is placed in different locations along the chain, is also observed. 

In our study of the semi-flexible ring chain we have observed a total lack of the n- or N- universality in g^ (f) . The 
correlation hole at small f is particularly sensitive on n even for chain separations n well exceeding the persistence 
length parameter cx A. Naturally, the correlation hole expands for small n due to the stiffness of the chain and 
consequently contracts for larger n. We may note also that in the latter case the behaviour of <?„ '(f) is somewhat 
similar to that of the chain with loose springs, leading to smaller values of 9. 

In the case of star polymers one has to distinguish the special role of the core monomer in addition to the end 
monomers, as well as the intra- and inter-arm correlations. Generally, in the coil state here we have a des Cloizeaux 
scaling behaviour with nearly universal 5, but with the exponent 9 being strongly sensitive on the monomer indices 
ij, number of arms /, as well as on the arm length (N — 1)//. A somewhat more stable behaviour for 9 was only seen 
for the intra-arm correlations when both monomers are fairly far along the arm from the core monomer, agreeing 
with the numerical values of #11=1,2 in an open chain. 

We have also studied the average shape of the above described polymers by means of computing the asphericity 
characteristics. We have found that the coil of an open chain is fairly aspherical and prolate (rod-like), while the coil 
of a ring is less so. Interestingly, the coil of a semi-stiff ring becomes more aspherical and prolate than that of the 
flexible ring as A increases or N decreases. Star polymers in the coil state are fairly aspherical and prolate for small 
numbers of arms /, becoming more spherical as / increases. Finally, the globule of a ring is nearly perfectly spherical. 

The static structure factor, related to the sum of the 3-d Fourier transforms of the correlation functions, has been 
investigated in terms of its appropriate scaling variables. Behaviour at small rescaled wave numbers p is universal, 
while the large p tail scales as a power law with the exponent related to the fractal dimension of the object. In 
cases of the flexible open chain, the ring coil and the globule the functions S(p) have shapes practically independent 
of polymer length N, which however is not true for semi-flexible polymers. Likewise, functions S(p) are strongly 
dependent on the number of arms / for star polymers. 

Finally, let us emphasise that we have only discussed the globule of a ring here because most of the results would 
be little, if at all, affected by the change of the connectivity matrix in the dense globular state, although this is a 
much more complicated issue meriting a separate study in case of fairy stiff chains, where, as we know, transitions 
into toroidal globules take placeEj. 

Thus, in this paper, we have seen that the degree of universality of various scaling functions and exponents is subject 
to more serious limitations than have been previously realised. The knowledge of various correlation functions of the 
polymer chain and their response to changes of various parameters provides us with the most detailed insight into 
the structure of the system. Development of theoretical methods, which could reproduce various features observed in 
this paper, is important for tackling larger system sizes and more complicated systems such as heteropolymers, for 
which simulations have severe limitations. This work is currently in progress and we hope to be able to return to such 
subject in the near future. 



ACKNOWLEDGMENTS 



The authors are grateful for interesting discussions to Professors F. Ganazzoli, T. Garel, J. Zinn- Justin, M. Daoud, 
B. Diinweg, and support from Enterprise Ireland: basic research grant SC/99/186 and international collaboration 
grants IC/2001/074 and BC/2001/034. 



1 J. des Cloizeaux, G. Jannink, Polymers in Solution (1990) Oxford Science Publ. 

2 H.H. Gan, B.C. Eu. J. Chem. Phys. 99 4084, 4103 (1993). 

3 M.P. Taylor, J.E.G. Lipson. J. Chem. Phys. 104, 4835 (1996); 106, 5181 (1997). 

4 A.M. Rubio, J.J. Freire, M. Bishop, J.H.R. Clarke. Macromol. 26 4018 (1993). 

5 J. P. Valleau. J. Chem. Phys. 104 3071 (1996). 

6 M. Bishop, J.H.R. Clarke. J. Chem. Phys. 94 3936 (1991). 
7 M. Bishop, J.H.R. Clarke. J. Chem. Phys. 95 4589 (1991). 
8 A. Baumgaitner. Z. Phys. B 42 265 (1981). 



10 



9 J. des Cloizeaux. J. Physique 41 223 (1980). 

10 B. Duplantier. J. Physique 48 569 (1986). 

11 B. Duplantier. J. Physique Lett. 46 L-751 (1985). 

12 Y. Oono, T. Ohta. Phys. Lett. 85A 480 (1981). 

13 L. Schafer, Excluded volume effects in polymer solutions as explained by the Renormalisation Group, Chapter 16, Springer- 
Verlag, Berlin (1999). 

14 R. Guida, J. Zinn-Justin, J. Phys. A 31, 8103 (1998). 

15 E.G. Timoshenko, Yu.A. Kuznetsov. Colloids and Surfaces A 190, 135 (2001). 

16 Yu.A. Kuznetsov, E.G. Timoshenko. J. Chem. Phys. Ill, 3744 (1999). 

17 F. Ganazzoli, Yu.A. Kuznetsov, E.G. Timoshenko. Macromol. Theory Simul. 10, 325 (2001). 

18 O. Kratky and G. Porod, Rec. Trav. Chim. 68, 1106 (1949); R.A. Harris and J.E. Hearst, J. Chem. Phys. 44, 2595 (1966). 

19 M.P. Allen and D.J. Tildesley (Ed.), Computer Simulation of Liquids (Clarendon Press, Oxford, 1987). 

20 K. Sole, W.H. Stockmayer. J. Chem. Phys. 54 2756 (1971). 

21 H.W. Diehl, E. Eisenriegler. J. Phys. A 22 L87 (1989). 

22 J. Rudnick, G. Gaspari. J. Phys. A 19 L191 (1986). 

23 O. Jagodzinski, E. Eisenriegler, K. Kremer. J. Phys. (France) I 2 2243 (1992). 

24 J.A. Aronovitz, D.R. Nelson, it J. Physique 47 1445 (1986). 

25 Yu.A. Kuznetsov, E.G. Timoshenko. II Nuovo Cimento 20D (12bis), 2265 (1998). 

26 J. -P. Hansen, I.R. McDonald, Theory of simple liquids (1990) Academic Press, London. 

27 J.C. Le Guillou, J. Zinn-Justin, Phys. Rev. Lett 39 95 (1977). 

28 E. Brezin, J.C. Le Guillou, J. Zinn-Justin, Phase Transitions and Critical Phenomena, Vol. 6, p. 125 (1976). 

29 W. H. Press, S. A. Teukolsky, W. T. Vetterling, and B. P. Flannery, Numerical Recipes in C (Cambridge University Press, 
1992). 

30 K. Schmitz, Dynamic light scattering by macromolecules, (1990) Academic Press, Boston. 



11 



FIGURE CAPTIONS 



FIG. 1. Plot of the mean-squared distances T> n (in units of £ 2 ) of flexible (A = 0) homopolymer rings in a good solvent 
(Vb = 0) vs the chain index n. Curves correspond to the values of the degree of polymerisation TV = 200, 300, 500 (from 
bottom to top). 

FIG. 2. Plot of the rescaled mean-squared distances Va — N~ 2v T> n (in units of £ 2 ) vs the rescaled chain index n — n/N 
of flexible homopolymer rings in a good solvent (Vb = 0, v = vw): bell-shaped curves for TV = 200 (solid line), TV = 300 
(long-dashed line), TV = 500 (short-dashed line); and in a poor solvent (Vb = 6fc_gT, v — 1/3): trapezoid-shaped curves for 
TV = 100 (solid line), TV = 150 (long-dashed line), TV = 200 (short-dashed line). 



FIG. 3. Plot of the rescaled correlation function g„ '(f) of homopolymer rings in a good solvent (Vb = 0) vs the rescaled 
spatial separation f corresponding to the chain index of half-ring (n = TV/2). Curves correspond to the values of the degree 
of polymerisation TV = 200 (solid line) and TV = 700 (pluses). Dashed line corresponds to half-ring n = TV/2 of TV = 200 
homopolymer, but with twice weaker connectivity constants «y = 1/2. 

FIG. 4. Plot of the mean-squared distances T> n (in units of I 2 ) of homopolymer rings in a poor solvent (Vb = 6 ksT) vs the 
chain index n. Curves correspond to the values of the degree of polymerisation TV = 50, 100, 150, 200 (from bottom to top). 

FIG. 5. Plot of the rescaled correlation function (f ) of homopolymer rings in a poor solvent (Vb = 6 ksT) vs the rescaled 
spatial separation f for the degree of polymerisation TV = 200. Curves correspond to the values of the chain index n = 10 
(short-dashed line), n = 15 (long-dashed line) and n = 100 (solid line). 

FIG. 6. Plot of the rescaled correlation function g^\f) of an open homopolymer in a good solvent (Vb = 0) vs the 
rescaled spatial separation r for the degree of polymerisation TV = 200. Curves correspond to the following values of the chain 
indices: end-end i,j = 1,200 (solid line and empty quadrangles), end-three quarters i,j = 1, 150 (long-dashed line and filled 
quadrangles), end-middle i,j — 1, 100 (short-dashed line and empty circles), and quarter-three quarters i,j = 50, 150 (dotted 
line and filled circles) from top to bottom in the left of the figure. Here lines have been obtained by fitting via Eq. mm and 
points correspond only to a subset of actual data to make the figure less cluttered. 



FIG. 7. Plot of the rescaled correlation function g„(r) of a semi-flexible (A/ft = 1) homopolymer ring in a good solvent 
(Vb = 0) vs the rescaled spatial separation f for the degree of polymerisation TV = 200. Curves correspond to the values of the 
chain index n = 10 (solid line), n = 20 (long-dashed line) and n = 50 (short -dashed line), and n = 100 (dotted line). Note 
also that these curves run from bottom to top in the left of the figure. 



FIG. 8. Plot of the rescaled correlation functions g? (r ) of a homopolymer star with the number of arms / = 12 and arm 
length (TV — 1)// = 50 in a good solvent (Vb = 0) vs the rescaled spatial separation f. Curves correspond to the following 
values of the chain indices: 0, a : m = 25 (solid line), 0, a : m = 50 (long-dashed line), a : n — 25, a : m — 50 (short-dashed 
line) and a : n = 50, b : m — 50 (dotted line). Note also that these curves run from bottom to top in the left of the figure. 



FIG. 9. Plot of the rescaled static structure factor S(p) of homopolymers vs the rescaled wave number p. Curves correspond 
to the following chains from top to bottom: open flexible coil with TV = 200 (long-dashed line), semi-flexible (A = 1) ring coil 
with TV = 200 (dotted line), semi-flexible (A = 1) ring coil with TV = 100 (short-dashed line), flexible ring coil with TV = 200 
(thick solid line), globule (Vb = 6ksT) of flexible ring with TV = 200 (solid circles curve), and the solid sphere with radius 
R s = y/(5/3)Hl (thin solid line). 
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FIG. 10. Plot of the rescaled static structure factor S(p) of star homopolymcrs vs the rescaled wave number p. Curves 
correspond to stars with the following numbers of arms from top to bottom: / = 3 (solid line), / = 6 (long-dashed line), / = 9 
(short-dashed line), / = 12 (dotted line). 
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TABLES 



TABLE I. Comparison of the exponents 8 and 9 between the results from Monte Carlo simulations and theoretical results 
(these are supplied with subscript theor) for ring and open homopolymer coils. Values in the second and fourth columns have 
been obtained by a four-parametric fit via Eq. (|l^) and the values in the fifth column (denoted as 9s fix have been obtained by a 
three-parametric fit via the same equation, but with fixed 8 = 8 t heor = 2.428363. The last column contains references to works 
with theoretical exponents. The notations in the first column follow des Cloizeaux convention: — end-end monomers, 1 — 
end-middle, 1' — end-three quarters, 1" — end-one quarter, and 2 — one quarter-three quarters of the chain respectively. 
Here and below reported errors are those from the fitting procedure only and do not necessarily account for statistical and 
other simulation errors. 





S 


$theor 


9 


8s fix 


Otheor 


Ref. 


Open TV = 200 


2.11 ±0.07 


2.428 ± 0.001 


0.36 ± 0.02 


0.30 ±0.02 


0.271 ±0.002 




1 


4 


Open TV = 200 1 


2.23 ±0.04 


2.428 ± 0.001 


0.56 ±0.01 


0.51 ±0.01 


~ 0.46 




1 




Open TV = 200 1' 


2.42 ±0.04 


2.428 ±0.001 


0.45 ±0.01 


0.45 ±0.01 


0.459 ±0.003 




1 




Open TV = 200 1" 


2.04 ±0.08 


2.428 ± 0.001 


0.68 ± 0.03 


0.56 ± 0.02 


~ 0.46 




1 




Open TV = 200 2 


2.39 ±0.07 


2.428 ± 0.001 


0.81 ±0.02 


0.80 ± 0.01 


0.71 ± 0.05 




1 




Ring TV = 300 


2.46 ±0.01 


2.428 ± 0.001 


0.79 ± 0.006 


0.816 ± 0.002 







TABLE II. Values of the exponents 8 and 9 for star polymers with / = 3,6, 12 arms and (TV — 1)// = 50 arm length in a 
good solvent. These have been obtained by a four-parametric fit via Eq. (hq) and the values for 9s fix have been obtained by 
a three-parametric fit via the same equation, but with fixed 8 = S t heor- The following notations for monomer pairs have been 
adopted: a : n,b : m, where a, b number arms and n, m number monomers within arms and refers to the core monomer. 







/ = 3 


/ = 6 


/ = 12 






8 


2.86 ±0.07 


2.66 ±0.04 


2.66 ±0.03 


0, a : m — 25 




8 


1.07 ±0.03 


2.01 ±0.03 


3.03 ±0.04 






9s fix 


1.25 ±0.02 


2.19 ±0.02 


3.30 ±0.02 






8 


2.54 ± 0.06 


2.57 ±0.04 


2.55 ±0.03 


0, a : m — 50 




9 


0.66 ±0.02 


1.05 ±0.02 


1.58 ±0.02 






9s fix 


0.70 ±0.01 


1.12 ±0.01 


1.66 ±0.01 






8 


2.31 ±0.06 


2.55 ±0.05 


2.39 ±0.03 


a : n = 25, a 


m = 50 


9 


0.57 ±0.02 


0.47 ±0.02 


0.54 ±0.01 






9s fix 


0.54 ±0.01 


0.51 ±0.01 


0.53 ±0.01 






8 


2.34 ± 0.06 


2.34 ± 0.03 


2.09 ±0.05 


a: n — 25, 6 : 


m = 25 


9 


0.84 ±0.03 


0.72 ±0.01 


0.59 ±0.03 






9s fix 


0.81 ±0.01 


0.69 ±0.01 


0.52 ±0.01 






8 


2.50 ±0.05 


2.47 ±0.03 


2.37 ±0.03 


a : n = 25, 6 : 


m = 50 


9 


0.41 ±0.01 


0.43 ±0.01 


0.36 ±0.01 






9s fix 


0.43 ±0.01 


0.44 ±0.01 


0.35 ±0.01 






8 


2.25 ±0.08 


2.43 ±0.04 


2.29 ±0.03 


a : n = 50, b : 


m = 50 


9 


0.35 ±0.03 


0.25 ±0.01 


0.22 ±0.01 






9s fix 


0.31 ±0.01 


0.26 ±0.01 


0.18 ±0.01 



TABLE III. Asphericity characteristics \ {a) (Tab. a), A3, S3, and A3, S3 (Tab. b) defined by Eqs. (|Io|[ll],|iJ|T|) respectively 
from Monte Carlo simulations. Here notations for chains correspond to: ring — flexible ring homopolymer coils, globule — 
flexible homopolymer rings in a poor solvent (Vb = 6ksT), open — flexible open homopolymer coils, stiff — semi-flexible 
ring homopolymer coils with stiffness constants A — 1,5, and star — flexible star homopolymers in a good solvent with arm 
length (TV - 1)// = 50. 
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a 



System 


A tlJ 


^ (2) 


\ { ' 


ring 








N = 100 


0.637 ± 0.001 


0.2637 ± 0.0008 


0.0995 ± 0.0004 


N = 200 


0.637 ± 0.001 


0.2643 ± 0.0008 


0.0986 ± 0.0004 


N = 300 


0.639 ± 0.001 


0.2651 ± 0.0008 


0.0957 ± 0.0004 


N = 500 


0.640 ± 0.001 


0.2654 ± 0.0008 


0.0945 ± 0.0004 


globule 








N = 100 


0.4230 ± 0.0004 


0.3242 ± 0.0003 


0.2529 ± 0.0001 


N = 150 


0.4098 ± 0.0002 


0.3265 ± 0.0002 


0.2637 ± 0.0001 


N = 200 


0.4021 ± 0.0002 


0.3277 ±0.0002 


0.2702 ± 0.0001 


open 








N = 150 


0.750 ±0.0015 


0.185 ±0.001 


0.0653 ± 0.0004 


N = 200 


0.7513 ± 0.0007 


0.1835 ±0.0005 


0.0652 ± 0.0002 


stiff A = 1 








N = 100 


0.6529 ± 0.0006 


0.2536 ± 0.0005 


0.0935 ± 0.0003 


N = 200 


0.6445 ± 0.0006 


0.2562 ± 0.0005 


0.0993 ± 0.0003 


N = 300 


0.6418 ± 0.0007 


0.2571 ± 0.0005 


0.1011 ±0.0003 


stiff A = 5 








N = 100 


0.6905 ± 0.0006 


0.2447 ± 0.0005 


0.0648 ± 0.0003 


at onn 

i\ — zuu 


U.Dt i-Z zt U.UUUD 


U.Z4 (O ± u.uuuo 


U.UolO =t u.uuuo 


JV = 300 


0.662 ± 0.001 


0.2537 ± 0.0005 


0.0841 ± 0.0003 


star 








/ = 3 


0.6752 ± 0.0009 


0.2451 ±0.0007 


0.0797 ± 0.0004 


/ = 6 


0.5492 ± 0.0007 


0.3036 ± 0.0005 


0.1472 ±0.0005 


/ = 9 


0.4990 ± 0.0006 


0.3163 ±0.0004 


0.1847 ±0.0004 


/ = 12 


0.4741 ± 0.0006 


0.3194 ±0.0003 


0.2065 ± 0.0004 
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b 



System 


A 3 


s 3 


As 


S 3 


ring 










N = 100 


0.253 ± 0.0015 


0.185 ±0.005 


0.282 ± 0.004 


0.243 ± 0.008 


N = 200 


0.254 ± 0.0016 


0.184 ± 0.005 


0.283 ± 0.004 


0.232 ± 0.008 


N = 300 


0.258 ± 0.0014 


0.187 ± 0.005 


0.286 ± 0.004 


0.236 ±0.008 


N = 500 


0.260 ±0.0015 


0.188 ± 0.005 


0.287 ± 0.004 


0.234 ± 0.008 


globule 










N = 100 


0.0255 ± 0.0001 


0.00317 ± 0.00007 


0.0261 ± 0.0002 


0.00344 ± 0.00009 


N = 150 


0.0192 ± 0.0001 


0.00173 ± 0.00004 


0.0188 ± 0.0002 


0.00184 ± 0.00007 


N = 200 


0.0155 ± 0.0001 


0.00118 ±0.00003 


0.0153 ± 0.0002 


0.00125 ± 0.00005 


open 










N = 150 


0.436 ± 0.0012 


0.550 ± 0.004 


0.548 ± 0.002 


0.906 ± 0.006 


N = 200 


0.4297 ± 0.0006 


0.539 ±0.002 


0.5492 ± 0.0008 


0.928 ±0.003 


stiff A = 1 










N = 100 


0.276 ± 0.001 


0.224 ± 0.002 


0.313 ± 0.002 


0.299 ± 0.004 


N = 200 


0.263 ± 0.001 


0.208 ± 0.002 


0.299 ± 0.002 


0.283 ± 0.004 


N = 300 


0.259 ± 0.001 


0.203 ± 0.002 


0.292 ± 0.002 


0.271 ±0.004 


stiff A = 5 










N = 100 


0.342 ± 0.001 


0.300 ± 0.002 


0.358 ± 0.002 


0.305 ± 0.004 


at onn 

i\ — zuu 


U.oU / zt U.UUl 


U.ZOO zt u.uuz 


U.ooD zt u.uuz 


n on/1 -i- n nn/i 
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N = 300 


0.293 ±0.001 


0.235 ± 0.002 


0.312 ±0.002 


0.219 ±0.004 


star 










/ = 3 


0.313 ±0.0014 


0.274 ± 0.003 


0.349 ± 0.003 


0.351 ±0.005 


/ = 6 


0.140 ±0.001 


0.0515 ±0.0007 


0.145 ±0.002 


0.057 ±0.001 


/ = 9 


0.0857 ± 0.0005 


0.0193 ±0.0005 


0.0874 ± 0.0008 


0.0197 ±0.0008 


/ = 12 


0.0622 ± 0.0003 


0.0112 ±0.0002 


0.0635 ± 0.0006 


0.0117 ±0.0003 
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